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Abstract 

Super Heavy quasi-stable particles are naturally produced in the early universe 
and could represent a substantial fraction of the Dark Matter: the so-called Super 
Heavy Dark Matter (SHDM). The decay of SHDM represents also a possible source 
of Ultra High Energy Cosmic Rays (UHECR), with a reliably calculated spectrum of 
the particles produced in the decay (oc E~ 19 ). The SHDM model for the production 
of UHECR can explain quantitatively only the excess of UHE events observed by 
AGASA. In the case of an observed spectrum not showing the AGASA excess the 
SHDM model can provide only a subdominant contribution to the UHECR flux. 
We discuss here the basic features of SHDM for the production of a subdominant 
UHECR flux, we refer our study to the possible signatures of the model at the 
Auger observatory discussing in particular the expected chemical composition and 
anisotropy. 

Key words: Super Heavy Dark Matter, Ultra High Energy Cosmic Rays, 
Anisotropy 



1 Introduction 



Ultra High Energy Cosmic Rays (UHECR) are the most energetic particles 
known in nature with energies exceeding 10 20 eV, the observation of particles 
with such high energies raises many interesting questions about their origin 
and composition that involve both astrophysics and particle physics. 
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Soon after the discovery of the Cosmic Microwave Background (CMB) radia- 
tion it was shown that the flux of UHE protons should be characterized by a 
sharp steepening at energy ~ 5 x 10 19 eV, due to the photo-pion production 
process on the CMB radiation field [T]. This effect is the well known Greisen- 
Zatsepin-Kuzmin (GZK) cut-off. After a few decades of observations the detec- 
tion of the GZK steepening of the CR flux is one of the major open problems 
in CR physics with experimental data still not conclusive. The AGASA ex- 
periment [2] observed 11 events with energy E > 10 20 eV in contrast with 
the expected depletion due to the GZK effect; on the other hand the HiRes 
experiment [3J observed an UHECR flux in agreement with the GZK cutoff 
[I]. The discrepancy between these two experimental results has been widely 
discussed and recently it has been shown that the GZK feature cannot be 
accurately determined with the small sample of events collected by AGASA 
and HiRes and that the discrepancy between the two experiments has a low 
statistical significance (at a 2a level) [5]. Moreover, AGASA and HiRes are 
based on different experimental techniques: a ground array the former and a 
fluorescence detector the latter. Taking this difference into account, the dif- 
fering results of AGASA and HiRes can be interpreted assuming the presence 
of a 30% systematic error in the relative energy determination. As shown by 
[dTIS] . correcting for these systematics the observations of AGASA and HiRes 
are brought into agreement at energies E < 10 20 eV and the discrepancy at 
the highest energies is softened. Recently, the Pierre Auger Observatory [TJ, 
under completion in Argentina, published its first observations of the UHECR 
flux combining the ground based and fluorescence detection techniques. This 
preliminary result seems to confirm the presence of the GZK steepening [S]. 

The presence of an excess of very high energy events, as claimed by AGASA, in- 
spired the introduction of several exotic models for the production of UHECR. 
These models, collectively called top-down, reproduce the excess of AGASA 
and give also an explanation for the lack of any clear astrophysical counterpart 
to the highest energy events observed. Indeed, at energies E > 5 x 10 19 eV 
the proton attenuation lenght is only about 20 4- 30 Mpc and an astrophysical 
source at this distance should have been seen at least in different frequency 
ranges. This evidence could imply that particles with energy E ~ 10 20 eV have 
a different (top-down) origin respect to those at lower energies. Many different 
ideas have been proposed among top-down models: strongly interacting neu- 
trinos j9] and new light hadrons [10] as unabsorbed signal carriers, Z-bursts 
|llj . Lorentz-invariance violation |12] , Topological Defects (TD) (see [13111] 
for a review), and Superheavy Dark Matter (SHDM) (see [15] for a review). 

The two models based on SHDM and TD have common features: in both 
cases UHECR are produced in the decay (SHDM) or annihilation (TD) of 
super-heavy particles, with a typical mass Mx > 10 13 GeV, that we will 
call collectively X-particles. As already discussed in [16], TD are distributed 
over cosmological distances therefore give only a marginal contribution to the 
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UHECR flux. We will not discuss this case here, concentrating our attention 
on the SHDM hypothesis. The possible existence of super heavy relic particles 
is an interesting conclusion of modern cosmology, being first suggested in 
connection with UHECR production [TTfTS] and later developed as a suitable 
candidate for Dark Matter (DM). 

Two main problems should be addressed in the discussion of SHDM models: 
how particles with very high mass (Mx > 10 13 GeV) can be quasi-stable, 
with a lifetime much longer than the age of the universe to, and how their 
abundance can be dominant in the universe today. The stability of SHDM 
can be achieved assuming the existence of a discrete gauge simmetry that 
protects the particle from decaying, in the same way as neutralino stability 
through R-parity in Super Symmetry (SUSY). This discrete simmetry can be 
weakly broken, assuring a lifetime Tx > to, through wormhole [17J or instanton 
[TS] effects, an example of a particle with a lifetime exceeding the age of the 
universe can be found in [19] . The abundance of SHDM can easily be dominant 
in the universe today, with a SHDM density Qshdm ~ ^dm- This effect can be 
obtained by gravitational production that resembles the production of density 
fluctuations during inflation [20] . 

The top-down hypothesis and, in particular, the UHECR production through 
SHDM decay can account only for the highest energy part of the observed 
spectrum, as in the case of AGASA excess at energies E > 5 x 10 19 eV [T6] . 
However, UHECR observations cannot exclude SHDM as explanation of the 
DM problem, assuming that SHDM is gravitationally produced than the X 
particle mass and density are unambiguously fixed, the only free parameter left 
to fit UHECR observations is the life-time tx- The observation of the AGASA 
excess fixes this value as tx — 10 20 y [T6|2T] . From the HiRes, Yakutsk and 
Auger data, that are compatible with the GZK steepening, follows only a 
lower bound on the value of t x . Within this lower bound it is still possible to 
test the SHDM hypothesis being connected with a subdominant contribution 
to the observed UHECR flux. In the present paper we will discuss such a 
possibility referring in particular to the observations of chemical composition 
and anisotropy of the Auger observatory. 

The UHECR production through the decay of SHDM shows three basic sig- 
natures that can be used to test the model: 

• SHDM particles, as any DM candidate, are clustered by gravitational inter- 
action and accumulated in the halo of our galaxy with a typical over-density 
of 5 ~ 2 x 10 5 [TTJ . Hence the UHECR spectrum from SHDM can overcome 
the constrain of the GZK steepening [17] . 

• In the decay of X-particles pions are extensively produced [21], therefore 
UHE neutrinos and photons are the dominant component of the primary 
flux [25. 
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• The non-central position of the Sun in the galactic halo results in an anisotropic 
flux of UHECR observed on earth [2TH22] . 

The photons (and neutrinos) dominated flux expected from SHDM is one 
of the most striking predictions of the model. Recently, the expected SHDM 
photon flux at energies E > 10 19 eV was compared with the photon fraction in 
the AGASA events [21J: this fraction was determined by different groups with 
different results, being around 50% in [23] and higher 67% in [21]. Another 
important analysis on the photon content in UHECR events was recently 
published by the Auger collaboration [25j. From this analysis, performed at 
the highest (E > 10 19 eV) energies, follows a stringent limit on the photon 
fraction that is around 2% at 10 19 eV and below 40% at 10 20 eV g5j. As 
recently discussed in [26], these limits do not exclude, but disfavor, the SHDM 
hypothesis for the explanation of the observed UHECR spectrum, in the next 
section we will come back to the expected photon fraction in SHDM models. 

The second most important signature of UHECR production by SHDM is the 
peculiar anisotropy expected. This anisotropy in the flux is guaranteed by 
our position respect to the Galactic Center (GC). In particular, the distance 
between the Sun and the outer boundary of the galaxy is larger in the GC 
direction respect to the anti-center direction. In order to evaluate the expected 
anisotropy it is necessary to assume a particular distribution of SHDM in the 
galaxy: in this respect numerical simulations of the galactic DM distribution 
show an increase of the DM density towards the GC as oc r _1 [27J or oc r~ L5 
[2"8] , this further enhancing the expected anisotropy. 

Already several authors have considered anisotropy computations, with re- 
liable predictions of the expected signal in the highest energy part of the 
spectrum [T3l21[22j . where the SHDM contribution can be dominant. The 
expected anisotropy was compared with the existing UHECR data taken in 
the northern hemisphere, this comparison reveals no contradiction between 
data and the predictions of the model [29j,[30j. On the other hand, the detec- 
tors located in the southern hemisphere are able to observe the GC and are 
much more sensitive to the expected anisotropy: in this respect the data of 
the old Sugar detector, in Australia, are only marginally consistent with the 
predictions of the SHDM model [3~T] . 

The Pierre Auger Observatory, due to its location in the southern emisphere 
and its large exposure, has a big potential for the detection of UHECR from 
the GC and the related anisotropy. In the present paper we will discuss such 
a potential, with particular attention to the lower energy region accessible to 
the detector, namely the range 3 x 10 18 4- 3 x 10 19 eV. At these energies the 
SHDM contribution to the flux, being subdominant, can be extracted from the 
data only using its peculiar anisotropy pattern. Although subdominant, UHE 
particles coming from the decay of SHDM are mainly photons that do not 
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suffer any isotropization process by magnetic fields, while 3 x 10 18 4- 3 x 10 19 
eV particles produced by astrophysical sources, being charged, are deflected 
by the galactic magnetic field [32J showing a substantially isotropic flux. In 
this sense, the low energy data of the Auger detector, with high statistics, can 
be used to shad light on the SHDM hypothesis, being any observed anisotropy 
directly connected only with SHDM. 

In order to compute the extragalactic UHECR flux we have used the Dip 
model, recently suggested to explain the behavior of UHECR at energies 
E > 10 18 eV ©. In this model the observed UHECR flux in the energy range 
10 18 eV 4 8 x 10 19 eV is explained as the effect of the combination of the adi- 
abatic and pair-production energy losses suffered by protons interacting with 
the CMB field. This interpretation of the observed spectra has important 
consequences mainly from the point of view of the origin and chemical com- 
position of UHECR, implying that, already at energies E > 10 18 eV, UHECR 
have an extragalactic origin and are mainly constituted by protons, with a 
steep injection spectrum oc E~~ f 7 > 2.5. 

The paper is organized as follows: in section [2] we will introduce the compu- 
tation of the expected flux from SHDM combining it with the extragalactic 
proton flux as in the Dip model, discussing in particular the expected chem- 
ical composition at the Auger observatory. In section [3] we will calculate the 
anisotropy of UHECR, referring our computation to the case of the Auger 
detector. Finally in section @] we will conclude discussing our results. 



2 Flux and chemical composition 

The first step to determine the UHECR flux produced in the decay of X par- 
ticles is the computation of the particles spectra. As it was already discussed 
in [12], this evaluation is particularly important because it represents a di- 
rect signature of the production mechanism that, in principle, can be detected 
experimentally. As discussed in the introduction, the mass of the decaying 
particle, Mx, that represents the total CMS energy y/s in the decay, is in the 
range 10 13 4 10 16 GeV. From the point of view of elementary particle physics 
the X-particle decay proceeds in a way similar to the e + e~ annihilation into 
hadrons, two or more off mass shell partons are produced initiating QCD 



2 The observed UHECR spectrum has also an alternative interpretation based on 
the so-called mixed composition scenario [33 . In this model UHECR with energies 
E > 10 19 eV are constituted by a mixture of extragalactic protons and heavy 
nuclei, with a flat injection spectrum 7 < 2.3, while, at lower energies E < 10 19 eV 
UHECR have still a galactic origin [33] . From the point of view of the SHDM induced 
anisotropy the "Dip" and "mixed composition" scenarios show similar features. 
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cascades. 



The determination of the particles energy spectra produced in the X-particle 
decay requires the extension of the existing QCD calculations for parton cas- 
cade from the TeV scale up to the Mx energy scale. The first calculations 
proposed in literature used an analytical approach based on the limiting spec- 
trum approximation [35] or an extension of the HERWIG [36] QCD Monte 
Carlo (MC) [3T] . More recently several papers appeared discussing the com- 
putation scheme involved in the SHDM decay. These schemes are based on two 
different approaches: the DGLAP evolution equations [T6l38|39j and a SUSY- 
QCD MC approach generalized to the high energies involved in the decay 
[KJffTB] . These computations predict quite accurately the secondary spectra 
from the X-particle decay and agree well with each other. Their most impor- 
tant outcome, concerning UHECR physics, is a flat energy spectrum of the 
produced particles with a behavior that can be approximated as dE/E 1,9 |16j . 

As already discussed in the introduction, an important signature of the SHDM 
model is represented by the peculiar UHECR composition, with a photon dom- 
inated flux. The exact determination of the photon fraction in the expected 
spectra is another outcome of the computation schemes discussed above, the 
typical photon/nucleon fraction expected in the cascade is r^/ N ~ 2 4- 3 [TB] . 
being only weekly dependent on the energy. 

The UHECR emissivity produced by the decays of SHDM in the halo of the 
galaxy can be simply written as: 

Ishdm(E) = —^—D(E)n x (R) (1) 



where D(E) is the energy spectrum (fragmentation function) of particles (pho- 
tons and protons GO) produced by the SHDM decay, and n x (R) is the density 
of X-particles in the galaxy as function of the distance R from the GC. In 
the present work we will assume that X-particles contribute to a substantial 
fraction of the galactic DM: therefore we will use the DM density as obtained 
by several numerical simulations. In particular we will use the Navarro- Frank- 
White (NFW) [27] and Moore [28] DM density profiles 

nx{R) = (R/R s nil R /R s r-* (2) 



with a = 1, 3/2 for the NFW and Moore case respectively, and R s = 45 Kpc 
as obtained in |29|. 



3 In the present paper we will not discuss neutrinos from SHDM, that is the subject 
of a forthcoming paper |41j 
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Fig. 1. The calculated spectrum of UHECRs from SHDM (dotted curve) in com- 
parison with the AGASA data (M x = 10 13 GeV left panel, M x = 10 14 GeV right 
panel). The spectrum from SHDM decays can explain only the highest energy events 
(AGASA excess). The continuos curve gives the spectrum of extragalactic protons 
(assuming the Dip model [6]). The sum of these two spectra is shown by thick 
continuos curve. 

Once determined the UHECR emissivity the expected flux can be simply 
computed as the line of sight integral of the emissivity ([1]) . Using the spherical 
symmetry of the DM density profile n x (R) we can identify the single line of 
sight simply by one angular variable that represents the angle formed by 
the axis Sun-GC and the line of sight itself. Therefore, the UHECR flux from 
the SHDM decay will be 

l^rnax ($) 

Jshdm{E,$) = —± D(E) I drn x {R{r)) (3) 

47T MvTx J 



where r is the line of sight coordinate and r max {d) is the distance to the 
boundary of the galactic halo in the d direction 

= R Q cos(tf) + ^R\ + R%sinH (4) 



where Rh = 100 Kpc is the radius of the galactic halo and R & = 8.5 Kpc is 
the distance Sun-GC. 



Let us now introduce the extragalactic (EG) component of the UHECR flux. 
As discussed in the introduction, we will use the Dip model [6] assuming an 
homogeneous distribution of sources and a proton dominated injection at the 
source. Under these hypothesis the UHECR flux can be computed as 



Jeg{E) 



c 

47T 



fdz 


dt 





dz 



Q inj (E g (E,z)) 



dE g (E,z) 
dE 



(5) 
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where: dt/dz = (Hq(1 + z)yft m (l + z) 3 + Q\) 1 represents the standard cos- 
mology (Q m = 0.27 and Q\ = 0.73 [12]); Qinj{E g ) oc _E~ 7 is the number of 
particles injected at the source per unit volume and energy, being the best fit 
value in the Dip model ^ g = 2.7 [6]; the generation energy E g (E,z) is deter- 
mined solving the losses equation dE g /dt = b(E g ,t) with the initial condition 
E g (t = t ) = E [6]; the quantity dE g /dE can be computed following the recipe 
given in [6lP£3] . 

The total UHECR flux will be the sum of the two contributions: galactic (from 
SHDM) and extragalactic (from homogeneous distributed sources) 

Juhecr{E, #) = Jeg(E) + J S hdm{E, #) (6) 



As it follows from equation ([3]) the flux from SHDM depends on two basic 
parameters: the X-particle life-time tx and mass Mx- While the value of Mx 
is directly connected with the production mechanism of SHDM, its life-time 
can be regarded as a free parameter that should be fixed in order to fit the 
observed UHECR events. The most remarkable creation mechanism for SHDM 
is gravitational production through time variable gravitational fields during 
the inflationary phase of the universe [H] . In this case the mass Mx fixes the 
present abundance of X-particles flsHDM- Assuming that SHDM constitute 
a substantial fraction of the DM, one has flsHDM — ^dm and to provide 
^shdm = 0.27, according to the WMAP observations [12], the mass of the 
X-particles must be in the range 10 13 4- 10 14 GeV. Fixing the X-particle mass 
M x = 10 13 , 10 14 GeV, in figured we plot the total UHECR flux © integrated 
over the AGASA field of view with the proper AGASA exposure. Fitting the 
AG AS A excess we obtain the life time value tx — 10 20 y The determination 
of the X-particle life time depends on the assumptions about the DM density 
profile with lower life-time in the case of the more concentrated Moore density 
profile (figured] refers to the NFW density profile). In figure[2]we have repeated 
the same computation of figured] but in the case of the Auger detector. In this 
case the observed flux [S] is almost consistent with the GZK suppression, hence 
the SHDM contribution to the flux is subdominant with a typical life-time of 
the order of Tx — 10 21 y. 

We shall now discuss the expected UHECR chemical composition in the frame- 
work of SHDM production. In figure [3] we have plotted the fraction of photons 
in the total flux computed under the same assumptions of figure [2] (Auger 
observations, in the two cases M x = 10 13 GeV, left panel, and M x = 10 14 
GeV, right panel). The result obtained, with a 2 4 3% photon content at 10 19 
is consistent with the Auger result that fixes the photon limit at 10 19 eV to 
about 2% [25|26] . the same result holds at 10 20 eV with an observed photon 
limit of 40% in the Auger data and an expected photon fraction, as follows 
from figure [3] of the same order. As in the case of the flux computation from 
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Fig. 2. The calculated spectrum of UHECRs from SHDM (dotted curve) in com- 
parison with the Auger data (M x = 10 13 GeV left panel, M x = 10 14 GeV right 
panel). The continuos curve gives the spectrum of extragalactic protons (assuming 
the Dip model [6j). The sum of these two spectra is shown by thick continuos curve. 

the chemical composition analysis, we can conclude that, although disfavored, 
the SHDM hypothesis is not ruled out by the available observations. 

Let us conclude this section discussing the extragalactic contribution to the 
flux due to SHDM. The dominant photon component of this flux is absorbed 
by the intergalactic medium, it is completely converted in electromagnetic cas- 
cades in a few Mpc, see the discussion in [13]. At the same time, the charged 
component (protons) suffers the interaction with CMB giving rise, at the high- 
est energies, to the photo-pion production process. These interactions deplete 
the high energy flux of the charged extragalactic SHDM contribution, making 
such flux vanishingly small if compared with the corresponding astrophysi- 
cal and galactic fluxes. For this reason we have neglected any extragalactic 
contribution to the SHDM flux. 



3 Anisotropy 



The anisotropy in the direction of the GC is another important prediction 
of the SHDM model for the production of UHECR. In particular, being the 
SHDM distributed in our galactic halo with a spherical symmetry around the 
GC, a pure dipole deviation from isotropy might be expected with the dipole 
vector pointing toward the GC itself. 

Independently of the top-down scenario, an anisotropy in the arrival directions 
of UHECR is also expected due to the cosmological Compton-Getting effect 
[injHT]. At energies E > 10 18 eV the UHECR flux is dominated by sources at 
cosmological distances, the movement of the Galaxy with respect to the CMB 
induces an anisotropy at 0.6% level, in any case smaller than the signals we 
will discuss here. 
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Fig. 3. Fraction of galactic SHDM photons in the total UHECR flux. The fluxes are 
computed in the case of the Auger detector assuming an NFW density profile and 
fixing M x = 10 13 GeV (left panel) and M x = 10 14 GeV (right panel). 

As already anticipated in the introduction, we are mainly interested in the 
anisotropy expected at the Auger observatory. To compute it we introduce 
the relative acceptance u of the detector, following [3B]- The number of the 
detected CR events are indeed distributed in the sky depending on both the 
real celestial anisotropy and the detector relative exposure ui. The quantity ui 
can be computed in terrestrial equatorial coordinates assuming that the detec- 
tor has a full-time operation, that is, a uniform exposure in right ascension a. 
For a detector at latitude A, fully efficient for particles arriving with a zenith 
angle 9 < 9 m (in the case of Auger, A ~ -35° and 9 m = 60° for E > 3 10 18 eV 
[38]). then it is possible to deduce the dependence of the relative acceptance 
on the declination 5: 



lu = ujq [cos(A) cos(<5) sin(a m ) + a m sin(A) sin(<5)] 



(7) 



where a m is given by 



a. 



if £ > 1 

7T iff<-l 
cos -1 (£) otherwise 



_ cos(6 m ) — sin(\)sin(5) 
cos(X)cos(8) 

and the constant ujq is determined by the integrated acceptance of the obser- 
vatory: 

2tt / u{6) cos(5)d5 = — - . 

J OA 

where Sa is the surface covered by the observatory, Aj = Sa^ sin 2 (# m ) is its 
integrated acceptance and (5 m in = o,Q — 9 m , 5 max = ao + 9 m ) is the observatory 
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field of view in declination. For the fully completed Auger detector, Sa = 3000 
km 2 and A T ~ 7000 km 2 sr for 6 m = 60°. 

Taking into account the Auger observatory acceptance function, the expected 
UHECR flux can then be computed as a function of the equatorial terrestrial 
coordinates (a, 5) 

Juhecr(E, a, 6) = J E g{E)uj(5) + Jshdm(E, a, S)u(S) (9) 

where, as in equation J EG is the extragalactic contribution [protons in the 
Dip model) and Jshdm is the galactic contribution from SHDMlj. The SHDM 
contribution to the flux is mainly composed by photons with a small admix- 
ture of protons (not larger than 30%). At the lowest energies E ~ 3 x 10 18 
eV protons are deflected by the galactic magnetic field and do not contribute 
to anisotropy. To take into account this effect and to bracket the possible 
anisotropy we have considered two extreme cases: (i) averaging, at all ener- 
gies, the SHDM proton contribution over the sky and (ii) not averaging this 
contribution. From the point of view of the computed anisotropy the accuracy 
of these two approaches depends on energy, the case (i) is most reliable at 
low energy (E < 3 x 10 19 eV) while case (ii) gives a better description at the 
highest energies (E > 3 x 10 19 eV). 

Using equation §§§ it is possible to determine the number of events at energy 
E > Eq expected to be collected during a time T . In order to apply an 
harmonic analysis, we integrate the events over 5, in 24 bins in a (i.e., Aa = 
15°): 

Oi+Aa S ma x 

N ai (> E ) = S A T J daJdE J J UHECR (E,a,S)cos(5)d5 (10) 

In figure H] and [5] we show, in the case of NF W and More respectively, the right 
ascension distribution of the expected events at the Auger observatory after 
T = 5 yr period of data taking: the four panels in each figure correspond to 
four different energies, E = 3 x 10 18 , 10 19 , 3 x 10 19 , 10 20 eV as labeled. In 
figure H] we have averaged the SHDM proton component over the sky while in 
figure [5] we kept the angular dependence of the SHDM proton flux. These two 
figures can be regarded as the two possible extreme cases of lowest (figure HJ) 
and highest anisotropy (figure E]) • 

To study the anisotropy we apply the Rayleigh's formalism (see e.g. [49] ) 

4 The relation between the angle defined in equation ([3|), and the galactic co- 
ordinates (b,l) is cos(i9) = cos(6)cos(£) and transforming to terrestrial equatorial 
coordinates (a, 5) one obtains Jshdm(E, a, 5) in equation ([9|). 
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Fig. 4. The expected UHECR events at the Auger observatory in 5 years data taking 
as function of right ascension. The SHDM contribution is obtained in the case of 
NFW density profile, assuming a SHDM mass of M x = 10 13 eV and averaging, at 
all energies, the SHDM proton component over the sky. 



to derive the amplitude and phase of the first harmonic in right ascension. 
Within this method the amplitude r and phase are given by: r = y/a 2 + b 2 
and (j) = atan(~), where a = -j| Y^iLi cos(aj), b = Y,i=x sm(cvj), and N=24 in 
our case. In tabled] we show the amplitude of the first harmonic in the case of 
the two density profiles and in the four different energy ranges. In table CD we 
have also quoted the two cases of minimum and maximum anisotropy obtained 
respectively averaging and not averaging the SHDM proton contribution over 
the sky. We do not list the phase because in all cases it always result to 
be = 17.7 hr (the galactic center coordinates being acc = 17h 42m 30s 
S GC = -28° 55' 18"). 

The quantity of interest is then the number of events, N evs , necessary to detect 
such amplitude and phase with a defined significance, S. Since the uncertainty 
on the amplitude is given by J j^—, we derive that N(S) = ^p. Thus to get 
a measurement of the expected anisotropies at 4 s.d. level, the number of 
required events is listed in the same table respectively for the two density 
profiles and in the case of minimum and maximum expected anisotropy. Such 
number of events roughly correspond to a maximum of 7 years and a minimum 
of 2 years of data taking in the Auger observatory. The Moore case indeed 
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Fig. 5. The expected UHECR events at the Auger observatory in 5 years data taking 
as function of right ascension. The SHDM contribution is obtained in the case of 
Moore density profile, assuming a SHDM mass of Mx = 10 13 eV and keeping, at 
all energies, the angular dependence of the SHDM proton component. 

corresponds to a more concentrated distribution of X-particles and thus gives 
a stronger signal from the GC, that will be detectable with less statistics. It is 
moreover important to stress that, in both density cases, already in the first 
energy bin under study, i.e., E Q = 3 x 10 18 eV, the anisotropy is at the level 
of few % and can be detectable in a few years by Auger. 

As already discussed the SHDM induced anisotropy is connected with the DM 
distribution in the galaxy, thus it is expected to appear on a typical scale that 
corresponds to the fiducial radius R s of this distribution as denned in equation 
(J2J). The fiducial radius R s as a typical value of 45 Kpc [29] and it corresponds 
to an expected anisotropy on a typical scale of 60° as follows from figures H] and 
|5j This is a large scale anisotropy if compared with the anisotropies already 
studied by Auger in the GC direction on a typical scale of a few degrees |50j . 

The result obtained is of particular importance because it shows the possibility 
of testing the SHDM hypothesis already at low energies where the Auger 
observatory has a larger statistic. Moreover, in the present evaluation we have 
bracketed the possible anisotropy taking into account the two cases of a SHDM 
proton component averaged and not averaged over the sky. 
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Energy (EeV) 


r (NFW-min) 


N(4 s.d.) 


r (Moore-min) 


N(4 s.d.) 


> 3 


0.015 


150000 


0.018 


100000 


> 10 


0.045 


16000 


0.060 


9000 


> 30 


0.100 


3200 


0.140 


1600 


> 100 


0.300 


350 


0.440 


150 


Energy (EeV) 


r (NFW-max) 


N(4 s.d.) 


r (Moore-max) 


N(4 s.d.) 


> 3 


0.019 


90000 


0.027 


45000 


> 10 


0.065 


7500 


0.092 


4000 


> 30 


0.150 


1400 


0.215 


700 


> 100 


0.480 


150 


0.690 


70 



Table 1 

Amplitude of first harmonic in right ascension, and number of events required to 
a detection at 4 standard deviations. Upper panel shows the case with minimum 
anisotropy (with isotropization of SHDM protons), lower panel shows the case of 
maximum anisotropy (with no isotropization of SHDM protons). 



The first case gives a reliable anisotropy at low energy (E < 10 19 eV) where 
the effect of the galactic magnetic field is important; the second case is more 
reliable at the highest energies (E >3x 10 19 eV) where the galactic magnetic 
field has no effect on the UHECR proton propagation. 



4 Conclusions 

In the present paper we have investigated the hypothesis of SHDM as a source 
of a subdominant component in the observed UHECR flux. The SHDM hy- 
pothesis is a viable candidate for the DM in the universe. Super heavy particles 
are naturally produced in time varying gravitational fields during the inflation 
era of the universe. In this context the observed DM density fixes the expected 
mass of SHDM in the range (10 13 GeV, 10 14 GeV). Moreover, SHDM can be 
quasi-stable, due to a super-weak discrete gauge symmetry breaking, with a 
life time exceeding the age of the universe. These characteristics of SHDM are 
of particular importance from the point of view of UHECR because they make 
SHDM a natural candidate for the production of ultra high energy particles. 
In order to explain the observed UHECR spectrum with SHDM the only free 
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parameter of the model is the SHDM particles life-time Tx- 

Super heavy particles are accumulated in the halo of our galaxy with a typical 
overdensity of 2.1 x 10 5 thus the corresponding UHECRs can be unabsorbed 
by CMB not showing the GZK steepening. Among UHECR observations only 
the AGASA data, with an excess of events violating the GZK suppression, 
can be used to fix a value of the SHDM life-time; from the HiRes, Yakutsk 
and Auger data, almost compatible with the GZK cut-off, follows only a lower 
bound of the SHDM life time. In any case these flux observations do not 
exclude the SHDM hypothesis and its role in the UHECR production can still 
be investigated as a subdominant contribution to the flux. 

The basic signatures of the SHDM model for the UHECR production are a 
flat injection spectrum oc E -1 ' 9 , a photon (and neutrino) dominated flux and 
an anisotropy due to the off-center position of the earth in the galaxy. In the 
present paper we have analyzed both the chemical composition and anisotropy 
expected at the Auger observatory, our basic findings can be summarized as 
follows. 

• The observed Auger spectrum is still compatible with a subdominant con- 
tribution from SHDM with a typical life-time of Tx > 10 21 y. 

• The chemical composition observed by Auger, with a photon limit around 
2% at 10 19 eV, is still compatible with the chemical composition expected 
in the SHDM model, that shows a photon content of the spectrum at the 
level of few % at 10 19 eV. 

• The anisotropy expected by SHDM can be detected by the Auger observa- 
tory and, even at the lowest energies, can be measured with a good accuracy. 

The anisotropy signal obtained in the present evaluation is of particular im- 
portance because it shows that already at low energy the SHDM hypothesis 
can be tested. At the lowest Auger energies E > 10 18 eV the charged com- 
ponent of the UHECR flux is isotropized by the galactic magnetic field while 
the photon dominated flux from SHDM does not. This is an important result 
connecting, at low energy, any observed anisotropy in the flux to the SHDM 
component. In the present paper we have performed a detailed computation 
of the expected anisotropy signal at Auger, showing how this signal is within 
reach of the Auger capabilities. 

Finally, at the highest energies, the anisotropy amplitude computed here 
should be regarded as an upper limit (see the discussion in section [3]). There- 
fore, the lack of detection by Auger in 5 or 7 years of full operation would not 
discard the SHDM hypothesis and, probably, a very large aperture new gen- 
eration experiment will be needed to settle down the question about UHECR 
produced by the decay of SHDM. 
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